Background

This analysis document compliments FIA NLS Models: Biomass vs. Stand Age. All of the background information from that document applies to these analyses, which are extensions to them. The difference between that document and this analysis is the use of different data subsets.

Here, we fit the models using: 1) a temporally-balanced dataset, where we take the first and most-recent plot record for all plots in the dataset, 2) a temporally-balanced dataset (same as #1), but which excludes plot locations which have experienced harvest (at any point over the study interval 2000-2022)

Below the model fitting procedure is implemented by ecoprovince:

Analysis 1: Temporally-balanced analysis

211 - Northeastern Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   4834     1203.6                                
## 2   4833     1092.6  1 111.05  491.22 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 52557.62
## 2     2 52091.41
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.45951    0.16992   2.704  0.00687 ** 
## alpha   0.84855    0.03503  24.221  < 2e-16 ***
## A     473.76484   37.00356  12.803  < 2e-16 ***
## k     205.69316   18.21698  11.291  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4755 on 4833 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 2.017e-06
##   (1 observation deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in Mod.Sel3 %in% c(1, "1a", "1b", "1c", 4) : 
##   object 'Mod.Sel3' not found
## Error in Mod.Sel3 %in% c(1, "1a", "1b", "1c", 4) : 
##   object 'Mod.Sel3' not found
## Error in Mod.Sel3 %in% c(1, "1a", "1b", "1c", 4) : 
##   object 'Mod.Sel3' not found
##   model      AIC
## 1     2 52091.41
## 2    2a       NA
## 3    2b       NA
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.45951    0.16992   2.704  0.00687 ** 
## alpha   0.84855    0.03503  24.221  < 2e-16 ***
## A     473.76484   37.00356  12.803  < 2e-16 ***
## k     205.69316   18.21698  11.291  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4755 on 4833 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 2.017e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   4832     1179.6                                
## 2   4831     1067.7  1 111.87  506.17 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 52091.41
## 2     4 52463.90
## 3     5 51983.93
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.58118    0.17752   3.274  0.00107 ** 
## alpha   0.84078    0.03411  24.646  < 2e-16 ***
## a      33.40070    1.99286  16.760  < 2e-16 ***
## b     118.88478    5.92661  20.059  < 2e-16 ***
## c     118.37512    5.85956  20.202  < 2e-16 ***
## d       1.01653    0.05233  19.425  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4701 on 4831 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.

plotting 2

212 - Laurentian Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1  12943     5165.4                                
## 2  12942     5044.1  1 121.33   311.3 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 135621.4
## 2     2 135315.7
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.95237    0.12759   7.464 8.91e-14 ***
## alpha   0.62169    0.03318  18.736  < 2e-16 ***
## A     165.96924    5.63575  29.449  < 2e-16 ***
## k      78.87707    3.15407  25.008  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6243 on 12942 degrees of freedom
## 
## Number of iterations to convergence: 4 
## Achieved convergence tolerance: 8.623e-06
##   (16 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1  12942     5044.1                                
## 2  12941     5021.8  1 22.266  57.378 3.839e-14 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 135315.7
## 2    2a 135260.4
## 3    2b 135317.6
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau   9.238e-01  1.259e-01   7.336 2.33e-13 ***
## alpha 6.214e-01  3.094e-02  20.082  < 2e-16 ***
## A     1.921e+02  9.346e+00  20.559  < 2e-16 ***
## k     1.085e+02  7.979e+00  13.593  < 2e-16 ***
## p     2.121e-02  2.492e-03   8.509  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6229 on 12941 degrees of freedom
## 
## Number of iterations to convergence: 9 
## Achieved convergence tolerance: 5.063e-06
##   (16 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1  12941     5088.9                                
## 2  12940     4912.4  1 176.51  464.95 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2a 135260.4
## 2     4 135432.3
## 3     5 134977.3
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.19051    0.13607   8.749   <2e-16 ***
## alpha   0.67333    0.02702  24.922   <2e-16 ***
## a      13.38590    0.51341  26.073   <2e-16 ***
## b      75.80154    2.36631  32.034   <2e-16 ***
## c     121.05709    5.45544  22.190   <2e-16 ***
## d       1.41595    0.04287  33.030   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6161 on 12940 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (16 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

221 - Eastern Broadleaf Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5442     948.80                                
## 2   5441     857.87  1 90.926  576.69 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 59550.63
## 2     2 59004.09
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.1644     0.1261   1.304    0.192    
## alpha   0.8050     0.0312  25.799   <2e-16 ***
## A     500.5614    30.8115  16.246   <2e-16 ***
## k     155.6108    11.1725  13.928   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3971 on 5441 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 2.272e-06
##   (1 observation deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value  Pr(>F)   
## 1   5441     857.87                             
## 2   5440     856.53  1 1.3429   8.529 0.00351 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 59004.09
## 2    2a 58997.56
## 3    2b 58973.93
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.23383    0.13009   1.797   0.0723 .  
## alpha   0.81624    0.03142  25.975   <2e-16 ***
## A     293.53650   21.44552  13.688   <2e-16 ***
## k      64.20469    6.41740  10.005   <2e-16 ***
## s       1.31952    0.05962  22.132   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3959 on 5440 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 2.275e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5440     940.37                                
## 2   5439     848.31  1 92.062  590.26 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 58973.93
## 2     4 59506.03
## 3     5 58947.03
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.24222    0.13018   1.861   0.0629 .  
## alpha   0.80822    0.03087  26.178  < 2e-16 ***
## a      21.08256    2.67854   7.871 4.22e-15 ***
## b     179.51639   10.45102  17.177  < 2e-16 ***
## c     156.83930   15.71326   9.981  < 2e-16 ***
## d       1.51879    0.08979  16.915  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3949 on 5439 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

222 - Midwest Broadleaf Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   3547     1178.6                                
## 2   3546     1136.5  1 42.068  131.25 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 38157.03
## 2     2 38030.00
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.03509    0.24350   4.251 2.18e-05 ***
## alpha   0.72686    0.05912  12.295  < 2e-16 ***
## A     445.56009   55.43529   8.037 1.24e-15 ***
## k     254.71684   34.82136   7.315 3.17e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5661 on 3546 degrees of freedom
## 
## Number of iterations to convergence: 8 
## Achieved convergence tolerance: 2.413e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value  Pr(>F)    
## 1   3546     1136.5                              
## 2   3545     1132.9  1 3.5914  11.238 0.00081 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 38030.00
## 2    2a 38020.77
## 3    2b       NA
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## 
## Parameters:
##         Estimate Std. Error t value Pr(>|t|)    
## tau     1.091545   0.248065   4.400 1.11e-05 ***
## alpha   0.753308   0.059393  12.683  < 2e-16 ***
## A     342.271688  43.589966   7.852 5.38e-15 ***
## k     170.045425  28.883552   5.887 4.29e-09 ***
## p      -0.015140   0.006514  -2.324   0.0202 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5653 on 3545 degrees of freedom
## 
## Number of iterations to convergence: 10 
## Achieved convergence tolerance: 5.957e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   3545     1160.1                                
## 2   3544     1108.8  1 51.282  163.91 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2a 38020.77
## 2     4 38104.80
## 3     5 37946.29
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.21145    0.25520   4.747 2.14e-06 ***
## alpha   0.76881    0.05439  14.134  < 2e-16 ***
## a      14.01898    1.39169  10.073  < 2e-16 ***
## b     101.92211    6.02465  16.918  < 2e-16 ***
## c     115.56990    8.21593  14.067  < 2e-16 ***
## d       1.19701    0.06667  17.955  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5593 on 3544 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

223 - Central Interior Broadleaf Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   6383     1208.6                                
## 2   6382     1143.1  1 65.521  365.81 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 66739.63
## 2     2 66385.70
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.80821    0.12030   6.718    2e-11 ***
## alpha   0.70093    0.03432  20.426   <2e-16 ***
## A     212.97304    8.46135  25.170   <2e-16 ***
## k      72.75583    3.94421  18.446   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4232 on 6382 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 2.25e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   6382     1143.1                                
## 2   6381     1142.5  1 0.6350  3.5469    0.0597 .  
## 3   6381     1136.2  0 0.0000                      
## 4   6380     1128.1  1 8.1539 46.1159 1.215e-11 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 66385.70
## 2    2a 66384.15
## 3    2b 66349.16
## 4    2c 66305.16
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.85759    0.12173   7.045 2.05e-12 ***
## alpha   0.71209    0.03381  21.062  < 2e-16 ***
## A     127.11413    4.61707  27.531  < 2e-16 ***
## k      37.02397    1.25609  29.475  < 2e-16 ***
## p       0.14094    0.01883   7.486 8.09e-14 ***
## s       2.17271    0.15279  14.220  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4205 on 6380 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 4.063e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   6381     1195.9                                
## 2   6380     1127.7  1 68.139   385.5 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2c 66305.16
## 2     4 66675.81
## 3     5 66303.16
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.84613    0.12109   6.988 3.08e-12 ***
## alpha   0.71122    0.03383  21.023  < 2e-16 ***
## a      17.80094    2.07031   8.598  < 2e-16 ***
## b      98.51432    3.90582  25.222  < 2e-16 ***
## c     113.18502    6.98395  16.206  < 2e-16 ***
## d       1.44570    0.07453  19.398  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4204 on 6380 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

231 - Southeastern Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7773     2698.5                                
## 2   7772     2568.4  1 130.03  393.48 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 85845.56
## 2     2 85463.52
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.24476    0.17019   7.314 2.86e-13 ***
## alpha   0.57095    0.02662  21.444  < 2e-16 ***
## A     257.12975   11.27279  22.810  < 2e-16 ***
## k      67.16853    3.01370  22.288  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5749 on 7772 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 1.572e-06
##   (14 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7772     2568.4                                
## 2   7771     2545.5  1 22.931  70.004 < 2.2e-16 ***
## 3   7771     2566.4  0  0.000                      
## 4   7770     2496.6  1 69.841 217.364 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 85463.52
## 2    2a 85395.79
## 3    2b 85459.40
## 4    2c 85246.85
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau   1.506e+00  1.836e-01   8.201 2.76e-16 ***
## alpha 6.512e-01  2.258e-02  28.840  < 2e-16 ***
## A     1.450e+02  6.248e+00  23.205  < 2e-16 ***
## k     3.016e+01  1.222e+00  24.686  < 2e-16 ***
## p     1.076e-01  7.674e-03  14.024  < 2e-16 ***
## s     1.963e+00  9.744e-02  20.143  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5668 on 7770 degrees of freedom
## 
## Number of iterations to convergence: 12 
## Achieved convergence tolerance: 5.033e-06
##   (14 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7771     2669.8                                
## 2   7770     2494.5  1 175.36  546.22 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2c 85246.85
## 2     4 85766.57
## 3     5 85240.29
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.47599    0.18195   8.112 5.75e-16 ***
## alpha   0.65085    0.02255  28.860  < 2e-16 ***
## a      16.12785    0.82771  19.485  < 2e-16 ***
## b     120.36292    5.58045  21.569  < 2e-16 ***
## c     112.83900    9.26790  12.175  < 2e-16 ***
## d       1.63350    0.06625  24.657  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5666 on 7770 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (14 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

232 - Outer Coastal Plain Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7905     4072.2                                
## 2   7904     3928.6  1 143.61  288.94 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 88878.43
## 2     2 88596.51
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.76811    0.18004   4.266 2.01e-05 ***
## alpha   0.55314    0.02987  18.521  < 2e-16 ***
## A     270.07655   14.26944  18.927  < 2e-16 ***
## k      71.69928    3.93397  18.226  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.705 on 7904 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 1.137e-06
##   (32 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7904     3928.6                                
## 2   7903     3859.0  1 69.597  142.53 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 88596.51
## 2    2a 88457.16
## 3    2b 88567.74
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau   9.272e-01  1.885e-01   4.920 8.85e-07 ***
## alpha 6.432e-01  2.456e-02  26.191  < 2e-16 ***
## A     4.083e+02  4.091e+01   9.981  < 2e-16 ***
## k     1.596e+02  2.131e+01   7.489 7.67e-14 ***
## p     2.506e-02  1.901e-03  13.185  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6988 on 7903 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 4.089e-06
##   (32 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7903     4029.5                                
## 2   7902     3773.8  1 255.69  535.39 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2a 88457.16
## 2     4 88799.14
## 3     5 88282.71
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.00852    0.19081   5.285 1.29e-07 ***
## alpha   0.71043    0.02138  33.233  < 2e-16 ***
## a      20.51866    1.02371  20.043  < 2e-16 ***
## b     121.27486    6.30094  19.247  < 2e-16 ***
## c     111.11452    9.25971  12.000  < 2e-16 ***
## d       1.50283    0.06832  21.998  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6911 on 7902 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (32 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

234 - Lower Mississippi Riverine Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    826     218.04                                
## 2    825     194.59  1 23.448  99.412 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 9094.121
## 2     2 9001.802
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau   3.865e-01  4.112e-01   0.940   0.3476    
## alpha 7.628e-01  6.984e-02  10.922   <2e-16 ***
## A     1.255e+03  5.318e+02   2.360   0.0185 *  
## k     4.893e+02  2.232e+02   2.192   0.0286 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4857 on 825 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 1.631e-06
##   (1 observation deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in nls(get(paste("f_", Mod.Sel1, "a", sep = "")), data = G_234,  : 
##   step factor 0.000488281 reduced below 'minFactor' of 0.000976562
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq Df Sum Sq F value  Pr(>F)  
## 1    825     194.59                            
## 2    823     192.48  2 2.1059  4.5022 0.01136 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 9001.802
## 2    2a       NA
## 3    2b       NA
## 4    2c 8996.781
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.39536    0.41113   0.962    0.337    
## alpha   0.78925    0.07021  11.242  < 2e-16 ***
## A     395.47217  189.29921   2.089    0.037 *  
## k      99.96670   70.69953   1.414    0.158    
## p      -0.01220    0.01423  -0.857    0.392    
## s       1.13109    0.24438   4.628 4.28e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4836 on 823 degrees of freedom
## 
## Number of iterations to convergence: 40 
## Achieved convergence tolerance: 3.1e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    824     216.69                                
## 2    823     192.29  1 24.399  104.43 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2c 8996.781
## 2     4 9092.966
## 3     5 8995.935
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.4024     0.4120   0.977  0.32908    
## alpha   0.7883     0.0701  11.246  < 2e-16 ***
## a       0.0000     5.9324   0.000  1.00000    
## b     340.8238   204.3728   1.668  0.09576 .  
## c     675.4071   940.5067   0.718  0.47288    
## d       2.5547     0.6884   3.711  0.00022 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4834 on 823 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

242 - Pacific Lowland Mixed Forest

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

251 - Prairie Parkland (Temperate)

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   1387     347.56                                
## 2   1386     336.63  1 10.934   45.02 2.834e-11 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 14526.30
## 2     2 14483.87
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.28977    0.25193   1.150     0.25    
## alpha   0.68210    0.09544   7.147 1.43e-12 ***
## A     261.45525   28.66195   9.122  < 2e-16 ***
## k     106.87482   14.79744   7.223 8.39e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4928 on 1386 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 9.101e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df     Sum Sq F value Pr(>F)
## 1   1386     336.63                             
## 2   1385     336.63  1 5.7789e-05   2e-04 0.9877
##   model      AIC
## 1     2 14483.87
## 2    2a 14485.87
## 3    2b 14469.11
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.34766    0.25540   1.361    0.174    
## alpha   0.70227    0.09501   7.391 2.51e-13 ***
## A     142.90894   13.48623  10.597  < 2e-16 ***
## k      39.02880    3.98802   9.787  < 2e-16 ***
## s       1.74387    0.18179   9.593  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.49 on 1385 degrees of freedom
## 
## Number of iterations to convergence: 8 
## Achieved convergence tolerance: 6.83e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   1385     337.23                                
## 2   1384     325.82  1 11.409  48.464 5.181e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 14469.11
## 2     4 14488.37
## 3     5 14442.53
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.34162    0.25201   1.356    0.175    
## alpha   0.70124    0.09431   7.435 1.82e-13 ***
## a      24.39435    4.05270   6.019 2.24e-09 ***
## b      93.28179    7.41989  12.572  < 2e-16 ***
## c     102.73710    8.92197  11.515  < 2e-16 ***
## d       1.13595    0.10865  10.455  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4852 on 1384 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

255 - Prairie Parkland (Subtropical)

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value  Pr(>F)  
## 1    439     237.09                            
## 2    438     235.09  1 2.0089  3.7429 0.05368 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 4609.958
## 2     2 4608.197
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)   
## tau    -0.1340     0.6019  -0.223  0.82393   
## alpha   0.3504     0.1744   2.009  0.04512 * 
## A     320.3752   106.7015   3.003  0.00283 **
## k     143.2259    52.7209   2.717  0.00686 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7326 on 438 degrees of freedom
## 
## Number of iterations to convergence: 8 
## Achieved convergence tolerance: 5.634e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1    438     235.09                         
## 2    437     234.74  1 0.3499  0.6514 0.4201
##   model      AIC
## 1     2 4608.197
## 2    2a 4609.539
## 3    2b       NA
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)   
## tau    -0.1340     0.6019  -0.223  0.82393   
## alpha   0.3504     0.1744   2.009  0.04512 * 
## A     320.3752   106.7015   3.003  0.00283 **
## k     143.2259    52.7209   2.717  0.00686 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7326 on 438 degrees of freedom
## 
## Number of iterations to convergence: 8 
## Achieved convergence tolerance: 5.634e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value   Pr(>F)   
## 1    437     226.64                              
## 2    436     221.69  1 4.9467  9.7286 0.001934 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 4608.197
## 2     4 4594.028
## 3     5 4586.274
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau    -0.2440     0.5501  -0.443 0.657650    
## alpha   0.5283     0.1577   3.350 0.000877 ***
## a       9.5645     1.9483   4.909 1.30e-06 ***
## b      88.5138    12.4743   7.096 5.24e-12 ***
## c      55.9913     5.9220   9.455  < 2e-16 ***
## d       1.0392     0.1138   9.129  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7131 on 436 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

261 - California Coastal Chaparral Forest and Shrub

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit

  • add s model: does not fit

  • add s+p model: does not fit

  • unable to fit model (only 64 observations)

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

262 - California Dry Steppe

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit

  • add s model: does not fit

  • add s+p model: does not fit

  • unable to fit model (0 observations)

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

263 - California Coastal Steppe - Mixed Forest and Redwood Forest

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

313 - Colorado Plateau Semi-Desert

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

315 - Southwest Plateau and Plains Dry Steppe and Shrub

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

321 - Chihuahuan Semi-Desert

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

322 - American Semidesert and Desert

model selection 1

## Error in if (any(nEQ <- vNms != make.names(vNms))) vNms[nEQ] <- paste0("`",  : 
##   missing value where TRUE/FALSE needed
## Error in if (any(nEQ <- vNms != make.names(vNms))) vNms[nEQ] <- paste0("`",  : 
##   missing value where TRUE/FALSE needed
##   model AIC
## 1     1  NA
## 2     2  NA
## Warning in min(AIC1_322$AIC, na.rm = T): no non-missing arguments to min;
## returning Inf
## Error in get(paste("nls_322.", Mod.Sel1, sep = "")) : 
##   object 'nls_322.' not found

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"
  • Cannot fit model
  • not enough data (only 3 observations)

331 - Great Plains/Palouse Dry Steppe

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"
  • Cannot fit model

332 - Great Plains Steppe

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value  Pr(>F)   
## 1    149     87.775                             
## 2    148     83.201  1 4.5741  8.1365 0.00496 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 1640.707
## 2     2 1634.572
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)   
## tau     3.0055     4.1286   0.728  0.46779   
## alpha   0.9714     0.2977   3.263  0.00137 **
## A     231.9254   226.9888   1.022  0.30857   
## k     219.8242   199.9830   1.099  0.27346   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7498 on 148 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 4.357e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in nls(get(paste("f_", Mod.Sel1, "a", sep = "")), data = G_332,  : 
##   number of iterations exceeded maximum of 50
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##   Res.Df Res.Sum Sq Df   Sum Sq F value Pr(>F)
## 1    148     83.201                           
## 2    147     83.108  1 0.093179  0.1648 0.6854
##   model      AIC
## 1     2 1634.572
## 2    2a       NA
## 3    2b 1636.402
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)   
## tau     3.0055     4.1286   0.728  0.46779   
## alpha   0.9714     0.2977   3.263  0.00137 **
## A     231.9254   226.9888   1.022  0.30857   
## k     219.8242   199.9830   1.099  0.27346   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7498 on 148 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 4.357e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: does not fit
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Error in nls(f_5, data = G_332, start = c(tau = tau.start, alpha = alpha.start,  : 
##   Convergence failure: iteration limit reached without convergence (10)
##   model      AIC
## 1     2 1634.572
## 2     4 1636.879
## 3     5       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)   
## tau     3.0055     4.1286   0.728  0.46779   
## alpha   0.9714     0.2977   3.263  0.00137 **
## A     231.9254   226.9888   1.022  0.30857   
## k     219.8242   199.9830   1.099  0.27346   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7498 on 148 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 4.357e-06
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: does not fit

predict and plot

plotting 2

341 - Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit

  • log-normal phi model: does not fit

  • model not fitted because only 62 observations

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

342 - Intermountain Semi-Desert

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

411 - Everglades

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

M211 - Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5103     963.25                                
## 2   5102     841.07  1 122.19  741.21 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 54262.24
## 2     2 53571.63
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.63454    0.16362   3.878 0.000107 ***
## alpha   0.80068    0.02686  29.811  < 2e-16 ***
## A     446.29921   29.09383  15.340  < 2e-16 ***
## k     198.11601   13.75314  14.405  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.406 on 5102 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 3.889e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5102     841.07                                
## 2   5101     833.07  1 8.0008   48.99 2.905e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 53571.63
## 2    2a 53524.82
## 3    2b       NA
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.51537    0.15576   3.309 0.000944 ***
## alpha   0.81328    0.02703  30.083  < 2e-16 ***
## A     324.77328   20.26162  16.029  < 2e-16 ***
## k     102.66625   10.74934   9.551  < 2e-16 ***
## p      -0.05130    0.01127  -4.551 5.47e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4041 on 5101 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 2.641e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5101     952.58                                
## 2   5100     828.03  1 124.55  767.12 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2a 53524.82
## 2     4 54209.36
## 3     5 53495.89
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.53028    0.15658   3.387 0.000713 ***
## alpha   0.81390    0.02671  30.470  < 2e-16 ***
## a      14.56503    2.98921   4.873 1.14e-06 ***
## b     155.65092   10.10200  15.408  < 2e-16 ***
## c     179.79825   19.79447   9.083  < 2e-16 ***
## d       1.58647    0.10396  15.261  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4029 on 5100 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M221 - Eastern Broadleaf Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5180     879.40                                
## 2   5179     824.61  1  54.79  344.11 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 57351.18
## 2     2 57019.77
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.69968    0.13285   5.267 1.45e-07 ***
## alpha   0.82371    0.04201  19.609  < 2e-16 ***
## A     266.46784   10.04611  26.524  < 2e-16 ***
## k      61.97132    3.14535  19.703  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.399 on 5179 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 3.081e-06
##   (3 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq Df  Sum Sq F value  Pr(>F)    
## 1   5179     824.61                               
## 2   5178     823.80  1  0.8171  5.1361 0.02347 *  
## 3   5178     817.30  0  0.0000                    
## 4   5177     805.02  1 12.2794 78.9674 < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 57019.77
## 2    2a 57016.63
## 3    2b 56975.59
## 4    2c 56899.12
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.87739    0.14139   6.206 5.87e-10 ***
## alpha   0.82469    0.03959  20.832  < 2e-16 ***
## A     162.45200    5.61416  28.936  < 2e-16 ***
## k      37.03352    1.05271  35.179  < 2e-16 ***
## p       0.20173    0.01784  11.310  < 2e-16 ***
## s       2.64792    0.19044  13.904  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3943 on 5177 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 1.479e-06
##   (3 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5178     862.65                                
## 2   5177     803.67  1 58.977  379.91 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2c 56899.12
## 2     4 57255.49
## 3     5 56890.45
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.87572    0.14097   6.212 5.64e-10 ***
## alpha   0.82430    0.03969  20.766  < 2e-16 ***
## a      31.76130    2.62008  12.122  < 2e-16 ***
## b     122.35089    4.61204  26.529  < 2e-16 ***
## c     103.12833    4.07801  25.289  < 2e-16 ***
## d       1.27367    0.05860  21.735  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.394 on 5177 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (3 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M223 - Ozark Broadleaf Forest Meadow

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    596     80.866                                
## 2    595     71.900  1 8.9663    74.2 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 6115.535
## 2     2 6047.140
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau    -0.07855    0.26298  -0.299    0.765    
## alpha   0.89957    0.09748   9.228  < 2e-16 ***
## A     298.68647   41.73938   7.156 2.45e-12 ***
## k      95.82372   18.87728   5.076 5.16e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3476 on 595 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 5.77e-06
##   (3 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in nls(get(paste("f_", Mod.Sel1, "c", sep = "")), data = G_M223,  : 
##   number of iterations exceeded maximum of 50
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df  Sum Sq F value Pr(>F)
## 1    595     71.900                          
## 2    594     71.708  1 0.19119  1.5837 0.2087
##   model      AIC
## 1     2 6047.140
## 2    2a 6047.545
## 3    2b 6048.497
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau    -0.07855    0.26298  -0.299    0.765    
## alpha   0.89957    0.09748   9.228  < 2e-16 ***
## A     298.68647   41.73938   7.156 2.45e-12 ***
## k      95.82372   18.87728   5.076 5.16e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3476 on 595 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 5.77e-06
##   (3 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    594     80.723                                
## 2    593     71.768  1 8.9544  73.988 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 6047.140
## 2     4 6118.473
## 3     5 6050.044
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau    -0.07855    0.26298  -0.299    0.765    
## alpha   0.89957    0.09748   9.228  < 2e-16 ***
## A     298.68647   41.73938   7.156 2.45e-12 ***
## k      95.82372   18.87728   5.076 5.16e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3476 on 595 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 5.77e-06
##   (3 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M231 - Ouachita Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    676     152.77                                
## 2    675     142.27  1 10.498  49.809 4.216e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 7025.553
## 2     2 6979.211
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.5697     0.5042   1.130    0.259    
## alpha   0.8583     0.1139   7.536 1.57e-13 ***
## A     315.8268    62.8438   5.026 6.44e-07 ***
## k     147.8427    33.1149   4.465 9.40e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4591 on 675 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 4.062e-06
##   (1 observation deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq Df    Sum Sq F value Pr(>F)
## 1    675     142.27                            
## 2    674     142.26  1 0.0098593  0.0467 0.8290
## 3    674     142.26  0 0.0000000               
## 4    673     142.26  1 0.0069224  0.0327 0.8564
##   model      AIC
## 1     2 6979.211
## 2    2a 6981.164
## 3    2b 6981.181
## 4    2c 6983.148
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.5697     0.5042   1.130    0.259    
## alpha   0.8583     0.1139   7.536 1.57e-13 ***
## A     315.8268    62.8438   5.026 6.44e-07 ***
## k     147.8427    33.1149   4.465 9.40e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4591 on 675 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 4.062e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Error in nls(f_5, data = G_M231, start = c(tau = tau.start, alpha = alpha.start,  : 
##   Convergence failure: iteration limit reached without convergence (10)
##   model      AIC
## 1     2 6979.211
## 2     4 7029.273
## 3     5       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.5697     0.5042   1.130    0.259    
## alpha   0.8583     0.1139   7.536 1.57e-13 ***
## A     315.8268    62.8438   5.026 6.44e-07 ***
## k     147.8427    33.1149   4.465 9.40e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4591 on 675 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 4.062e-06
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: does not fit

predict and plot

plotting 2

M242 - Cascade Mixed Forest

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

M261 - Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value  Pr(>F)  
## 1    326    100.337                            
## 2    325     98.647  1 1.6905  5.5696 0.01887 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 3504.259
## 2     2 3500.669
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     1.5690     1.1894   1.319  0.18805    
## alpha   0.6751     0.2661   2.537  0.01165 *  
## A      85.0151    18.2427   4.660 4.61e-06 ***
## k      22.9171     7.8626   2.915  0.00381 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5509 on 325 degrees of freedom
## 
## Number of iterations to convergence: 12 
## Achieved convergence tolerance: 7.882e-06
##   (1 observation deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq Df  Sum Sq F value Pr(>F)  
## 1    325     98.647                            
## 2    324     97.539  1 1.10741  3.6785 0.0560 .
## 3    324     93.712  0 0.00000                 
## 4    323     93.308  1 0.40411  1.3989 0.2378  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 3500.669
## 2    2a 3498.954
## 3    2b 3485.785
## 4    2c 3486.363
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     1.8945     1.2923   1.466   0.1436    
## alpha   0.6580     0.2763   2.381   0.0178 *  
## A      65.9248    13.3842   4.926 1.34e-06 ***
## k      30.8188     3.7464   8.226 4.75e-15 ***
## s      42.5340   187.7406   0.227   0.8209    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5378 on 324 degrees of freedom
## 
## Number of iterations to convergence: 21 
## Achieved convergence tolerance: 8.771e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value  Pr(>F)  
## 1    324     96.856                            
## 2    323     95.476  1   1.38  4.6688 0.03145 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 3485.785
## 2     4 3496.641
## 3     5 3493.920
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     1.8945     1.2923   1.466   0.1436    
## alpha   0.6580     0.2763   2.381   0.0178 *  
## A      65.9248    13.3842   4.926 1.34e-06 ***
## k      30.8188     3.7464   8.226 4.75e-15 ***
## s      42.5340   187.7406   0.227   0.8209    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5378 on 324 degrees of freedom
## 
## Number of iterations to convergence: 21 
## Achieved convergence tolerance: 8.771e-06
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M262 - California coastal range - coniferous forest - open woodland - shrub meadow

Model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit

  • log-normal phi model: does not fit

  • model can fit - but K is negative (only 19 observations) - model excluded

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M313 - Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M331 - Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M332 - Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M333 - Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M334 - Black Hills Coniferous Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    299     124.35                                
## 2    298     113.37  1 10.974  28.846 1.581e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 3013.807
## 2     2 2987.903
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau    0.07662    0.66159   0.116  0.90788    
## alpha  0.77829    0.12977   5.998 5.79e-09 ***
## A     99.54181   19.23214   5.176 4.18e-07 ***
## k     53.05797   17.65911   3.005  0.00289 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6168 on 298 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 6.942e-06
##   (4 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq Df   Sum Sq F value Pr(>F)
## 1    298     113.37                           
## 2    297     113.26  1 0.115630  0.3032 0.5823
## 3    297     113.00  0 0.000000               
## 4    296     112.91  1 0.094618  0.2480 0.6188
##   model      AIC
## 1     2 2987.903
## 2    2a 2989.595
## 3    2b 2988.916
## 4    2c 2990.663
## Warning in `[<-.data.frame`(`*tmp*`, nls.param.df_bal$Code == "M334", , :
## provided 33 variables to replace 32 variables
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau    0.07662    0.66159   0.116  0.90788    
## alpha  0.77829    0.12977   5.998 5.79e-09 ***
## A     99.54181   19.23214   5.176 4.18e-07 ***
## k     53.05797   17.65911   3.005  0.00289 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6168 on 298 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 6.942e-06
##   (4 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    297     124.07                                
## 2    296     113.29  1 10.775  28.151 2.203e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 2987.903
## 2     4 3017.126
## 3     5 2991.690
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau    0.07662    0.66159   0.116  0.90788    
## alpha  0.77829    0.12977   5.998 5.79e-09 ***
## A     99.54181   19.23214   5.176 4.18e-07 ***
## k     53.05797   17.65911   3.005  0.00289 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6168 on 298 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 6.942e-06
##   (4 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M341 - Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2


Fitted parameters

Best / selected models by ecoprovince

Code Ecoregion Sel.Mod.2 Sel.Mod.3 Best.Mod
211 Northeastern Mixed Forest 2 5 5
212 Laurentian Mixed Forest 2a 5 5
221 Eastern Broadleaf Forest 2b 5 5
222 Midwest Broadleaf Forest 2a 5 5
223 Central Interior Broadleaf Forest 2c 5 5
231 Southeastern Mixed Forest 2c 5 5
232 Outer Coastal Plain Mixed Forest 2a 5 5
234 Lower Mississippi Riverine Forest 2c 5 5
242 Pacific Lowland Mixed Forest NA NA NA
251 Prairie Parkland (Temperate) 2b 5 5
255 Prairie Parkland (Subtropical) 2 5 5
261 California Coastal Chaparral Forest and Shrub NA NA NA
262 California Dry Steppe NA NA NA
263 California Coastal Steppe - Mixed Forest and Redwood Forest NA NA NA
313 Colorado Plateau Semi-Desert NA NA NA
315 Southwest Plateau and Plains Dry Steppe and Shrub NA NA NA
321 Chihuahuan Semi-Desert NA NA NA
322 American Semidesert and Desert NA NA NA
331 Great Plains/Palouse Dry Steppe NA NA NA
332 Great Plains Steppe 2 2 2
341 Intermountain Semi-Desert and Desert NA NA NA
342 Intermountain Semi-Desert NA NA NA
411 Everglades NA NA NA
M211 Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow 2a 5 5
M221 Central Appalachian Broadleaf Forest - Coniferous Forest - Meadow 2c 5 5
M223 Ozark Broadleaf Forest Meadow 2 2 2
M231 Ouachita Mixed Forest 2 2 2
M242 Cascade Mixed Forest NA NA NA
M261 Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow 2b 2b 2b
M262 California Coastal Range Coniferous Forest - Open Woodland - Shrub - Meadow NA NA NA
M313 Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow NA NA NA
M331 Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow NA NA NA
M332 Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow NA NA NA
M333 Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow NA NA NA
M334 Black Hills Coniferous Forest 2 2 2
M341 Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow NA NA NA

table by ecoprovince

Code Ecoregion region n.obs n.plots tau tau.variance tau.2.5 tau.97.5 alpha alpha.variance alpha.2.5 alpha.97.5 A A.2.5 A.97.5 k k.2.5 k.97.5 a a.2.5 a.97.5 b b.se b.2.5 b.97.5 c c.2.5 c.97.5 d d.2.5 d.97.5
211 Northeastern Mixed Forest east 4838 2419 0.5811807 0.0315130 0.2331627 0.9291986 0.8407842 0.0011638 0.7739056 0.9076629 473.76484 401.22103 546.30865 205.69316 169.97959 241.40674 33.400705 29.493797 37.30761 118.88478 NA 107.26593 130.50363 118.37512 106.88772 129.86252 1.016528 0.9139327 1.119122
212 Laurentian Mixed Forest east 12962 6481 1.1905141 0.0185145 0.9238008 1.4572274 0.6733287 0.0007299 0.6203705 0.7262870 192.14309 173.82377 210.46242 108.46141 92.82129 124.10153 13.385905 12.379552 14.39226 75.80154 NA 71.16321 80.43987 121.05709 110.36363 131.75055 1.415947 1.3319180 1.499977
221 Eastern Broadleaf Forest east 5446 2723 0.2422230 0.0169479 -0.0129899 0.4974359 0.8082238 0.0009532 0.7476984 0.8687492 293.53650 251.49470 335.57830 64.20469 51.62403 76.78536 21.082556 15.831554 26.33356 179.51639 NA 159.02822 200.00457 156.83930 126.03502 187.64358 1.518786 1.3427586 1.694814
222 Midwest Broadleaf Forest east 3552 1776 1.2114537 0.0651252 0.7111070 1.7118003 0.7688102 0.0029588 0.6621620 0.8754583 342.27169 256.80774 427.73563 170.04542 113.41537 226.67548 14.018976 11.290384 16.74757 101.92211 NA 90.10999 113.73424 115.56990 99.46147 131.67834 1.197007 1.0662985 1.327715
223 Central Interior Broadleaf Forest east 6388 3194 0.8461273 0.0146625 0.6087525 1.0835021 0.7112156 0.0011445 0.6448959 0.7775352 127.11413 118.06312 136.16514 37.02397 34.56160 39.48633 17.800944 13.742432 21.85946 98.51432 NA 90.85761 106.17104 113.18502 99.49412 126.87591 1.445696 1.2995960 1.591796
231 Southeastern Mixed Forest east 7790 3895 1.4759858 0.0331065 1.1193112 1.8326605 0.6508476 0.0005086 0.6066399 0.6950552 144.97930 132.73182 157.22678 30.16077 27.76581 32.55574 16.127853 14.505323 17.75038 120.36292 NA 109.42373 131.30210 112.83900 94.67141 131.00659 1.633502 1.5036347 1.763370
232 Outer Coastal Plain Mixed Forest east 7940 3970 1.0085159 0.0364082 0.6344793 1.3825526 0.7104306 0.0004570 0.6685254 0.7523358 408.34768 328.15225 488.54310 159.56714 117.80214 201.33213 20.518658 18.511922 22.52539 121.27486 NA 108.92335 133.62637 111.11452 92.96303 129.26601 1.502835 1.3689187 1.636751
234 Lower Mississippi Riverine Forest east 830 415 0.4023597 0.1697654 -0.4063858 1.2111052 0.7883253 0.0049140 0.6507290 0.9259216 395.47217 23.90610 767.03825 99.96670 -38.80593 238.73932 0.000000 -11.644423 11.64442 340.82383 NA -60.32947 741.97714 675.40707 -1170.66705 2521.48119 2.554711 1.2035154 3.905907
242 Pacific Lowland Mixed Forest west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
251 Prairie Parkland (Temperate) east 1392 696 0.3416233 0.0635084 -0.1527370 0.8359836 0.7012428 0.0088949 0.5162318 0.8862538 142.90894 116.45330 169.36458 39.02880 31.20558 46.85201 24.394348 16.444251 32.34445 93.28179 NA 78.72633 107.83724 102.73710 85.23504 120.23915 1.135955 0.9228090 1.349101
255 Prairie Parkland (Subtropical) east 444 222 -0.2439613 0.3026388 -1.3251902 0.8372676 0.5283394 0.0248661 0.2184128 0.8382661 320.37518 110.66455 530.08580 143.22594 39.60855 246.84333 9.564518 5.735244 13.39379 88.51384 NA 63.99654 113.03113 55.99135 44.35208 67.63061 1.039193 0.8154724 1.262914
261 California Coastal Chaparral Forest and Shrub west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
262 California Dry Steppe west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
263 California Coastal Steppe - Mixed Forest and Redwood Forest west 4 2 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
313 Colorado Plateau Semi-Desert west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
315 Southwest Plateau and Plains Dry Steppe and Shrub west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
321 Chihuahuan Semi-Desert west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
322 American Semidesert and Desert west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
331 Great Plains/Palouse Dry Steppe west 118 59 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
332 Great Plains Steppe west 154 77 3.0054586 17.0455900 -5.1532211 11.1641382 0.9713883 0.0886175 0.3831227 1.5596540 231.92537 -216.63238 680.48311 219.82421 -175.36661 615.01504 NA NA NA NA NA NA NA NA NA NA NA NA NA
341 Intermountain Semi-Desert and Desert west 4 2 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
342 Intermountain Semi-Desert west 2 1 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
411 Everglades east 66 33 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M211 Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow east 5108 2554 0.5302798 0.0245169 0.2233184 0.8372412 0.8138994 0.0007135 0.7615330 0.8662658 324.77328 285.05180 364.49475 102.66625 81.59293 123.73957 14.565033 8.704894 20.42517 155.65092 NA 135.84666 175.45518 179.79825 140.99259 218.60392 1.586471 1.3826673 1.790274
M221 Central Appalachian Broadleaf Forest - Coniferous Forest - Meadow east 5186 2593 0.8757158 0.0198734 0.5993491 1.1520824 0.8243016 0.0015756 0.7464844 0.9021187 162.45200 151.44589 173.45812 37.03352 34.96976 39.09727 31.761299 26.624837 36.89776 122.35089 NA 113.30935 131.39243 103.12833 95.13370 111.12296 1.273667 1.1587869 1.388547
M223 Ozark Broadleaf Forest Meadow east 602 301 -0.0785477 0.0691598 -0.5950347 0.4379392 0.8995732 0.0095032 0.7081176 1.0910287 298.68647 216.71203 380.66091 95.82372 58.74951 132.89792 NA NA NA NA NA NA NA NA NA NA NA NA NA
M231 Ouachita Mixed Forest east 680 340 0.5697483 0.2542625 -0.4203280 1.5598245 0.8582543 0.0129720 0.6346241 1.0818845 315.82677 192.43391 439.21963 147.84274 82.82219 212.86329 NA NA NA NA NA NA NA NA NA NA NA NA NA
M242 Cascade Mixed Forest west 34 17 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M261 Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow west 330 165 1.8945480 1.6701111 -0.6478649 4.4369609 0.6580466 0.0763610 0.1144095 1.2016837 65.92482 39.59387 92.25577 30.81882 23.44844 38.18920 NA NA NA NA NA NA NA NA NA NA NA NA NA
M262 California Coastal Range Coniferous Forest - Open Woodland - Shrub - Meadow west 8 4 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M313 Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M331 Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M332 Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow west 20 10 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M333 Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow west 22 11 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M334 Black Hills Coniferous Forest west 306 153 0.0766154 0.4376976 -1.2253593 1.3785901 0.7782879 0.0168399 0.5229091 1.0336667 99.54181 61.69380 137.38982 53.05797 18.30561 87.81033 NA NA NA NA NA NA NA NA NA NA NA NA NA
M341 Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA

plot tau

map

## OGR data source with driver: ESRI Shapefile 
## Source: "C:\Users\hogan.jaaron\Dropbox\FIA_R\Mapping\S_USA.EcoMapProvinces\S_USA.EcoMapProvinces.shp", layer: "S_USA.EcoMapProvinces"
## with 37 features
## It has 17 fields
## Integer64 fields read as strings:  PROVINCE_ PROVINCE_I

plot alpha (biomass compensation effect)

plot A (asymptote of B)

## Warning: Removed 19 rows containing missing values (`geom_point()`).

plot k (stand age at half biomass asymptote)

## Warning: Removed 19 rows containing missing values (`geom_point()`).


Analysis 2: No-harvest

211 - Northeastern Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   4767     1100.1                                
## 2   4766     1043.0  1 57.113  260.99 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 51935.28
## 2     2 51682.97
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.55012    0.29122   5.323 1.07e-07 ***
## alpha   1.10900    0.06402  17.322  < 2e-16 ***
## A     358.60119   28.52883  12.570  < 2e-16 ***
## k     179.16216   14.20596  12.612  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4678 on 4766 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 2.058e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value   Pr(>F)   
## 1   4766       1043                              
## 2   4765       1041  1 1.9622   8.982 0.002741 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 51682.97
## 2    2a 51675.99
## 3    2b 51658.42
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.45638    0.28154   5.173 2.40e-07 ***
## alpha   1.12704    0.06463  17.440  < 2e-16 ***
## A     214.23259   19.48572  10.994  < 2e-16 ***
## k      72.38593    8.79876   8.227 2.46e-16 ***
## s       1.30993    0.06714  19.511  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4665 on 4765 degrees of freedom
## 
## Number of iterations to convergence: 8 
## Achieved convergence tolerance: 8.183e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   4765     1083.1                                
## 2   4764     1027.3  1 55.751  258.53 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 51658.42
## 2     4 51865.15
## 3     5 51615.07
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.58266    0.29313   5.399 7.02e-08 ***
## alpha   1.10888    0.06432  17.240  < 2e-16 ***
## a      23.88706    2.03910  11.714  < 2e-16 ***
## b     108.49559    7.02793  15.438  < 2e-16 ***
## c     132.60613    8.82813  15.021  < 2e-16 ***
## d       1.19405    0.06821  17.506  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4644 on 4764 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

212 - Laurentian Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1  17134     6198.8                                
## 2  17133     6050.3  1  148.5  420.51 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 178725.1
## 2     2 178311.6
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.10028    0.13197   8.337   <2e-16 ***
## alpha   0.81927    0.03764  21.768   <2e-16 ***
## A     152.07486    4.72755  32.168   <2e-16 ***
## k      73.19820    2.56503  28.537   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5943 on 17133 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 8.52e-06
##   (13 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value   Pr(>F)   
## 1  17133     6050.3                              
## 2  17132     6047.1  1 3.2808  9.2947 0.002302 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 178311.6
## 2    2a 178304.3
## 3    2b 178234.1
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.16176    0.13429   8.651   <2e-16 ***
## alpha   0.83619    0.03798  22.015   <2e-16 ***
## A     112.92651    4.28367  26.362   <2e-16 ***
## k      41.69107    2.16303  19.274   <2e-16 ***
## s       1.30455    0.03898  33.467   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5929 on 17132 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 9.046e-06
##   (13 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1  17132     6118.8                                
## 2  17131     5964.1  1  154.7  444.35 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 178234.1
## 2     4 178506.3
## 3     5 178069.5
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.21466    0.13596   8.934   <2e-16 ***
## alpha   0.82392    0.03678  22.402   <2e-16 ***
## a      13.98018    0.80587  17.348   <2e-16 ***
## b      71.24841    2.09216  34.055   <2e-16 ***
## c     116.92499    4.25601  27.473   <2e-16 ***
## d       1.41175    0.04295  32.869   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.59 on 17131 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (13 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

221 - Eastern Broadleaf Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5798     981.08                                
## 2   5797     926.82  1 54.253  339.34 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 63460.21
## 2     2 63132.20
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.51229    0.15326   3.343 0.000835 ***
## alpha   0.88535    0.04543  19.490  < 2e-16 ***
## A     508.99853   32.76656  15.534  < 2e-16 ***
## k     174.84347   12.40713  14.092  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3998 on 5797 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 5.849e-06
##   (1 observation deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5797     926.82                                
## 2   5796     919.97  1 6.8551  43.189 5.405e-11 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 63132.20
## 2    2a 63091.14
## 3    2b 58973.93
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.23383    0.13009   1.797   0.0723 .  
## alpha   0.81624    0.03142  25.975   <2e-16 ***
## A     293.53650   21.44552  13.688   <2e-16 ***
## k      64.20469    6.41740  10.005   <2e-16 ***
## s       1.31952    0.05962  22.132   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3959 on 5440 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 2.275e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   5796     965.51                                
## 2   5795     910.77  1 54.739  348.29 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 58973.93
## 2     4 63371.45
## 3     5 63034.87
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.23383    0.13009   1.797   0.0723 .  
## alpha   0.81624    0.03142  25.975   <2e-16 ***
## A     293.53650   21.44552  13.688   <2e-16 ***
## k      64.20469    6.41740  10.005   <2e-16 ***
## s       1.31952    0.05962  22.132   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3959 on 5440 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 2.275e-06
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

222 - Midwest Broadleaf Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   4359     1525.4                                
## 2   4358     1473.0  1 52.402  155.04 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 46875.15
## 2     2 46724.67
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.79200    0.23894   3.315 0.000925 ***
## alpha   0.84551    0.06349  13.318  < 2e-16 ***
## A     630.56526   94.45976   6.675 2.77e-11 ***
## k     370.30471   59.25802   6.249 4.53e-10 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5814 on 4358 degrees of freedom
## 
## Number of iterations to convergence: 8 
## Achieved convergence tolerance: 7.781e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   4358     1473.0                                
## 2   4357     1460.9  1  12.03  35.879 2.269e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 46724.67
## 2    2a 46690.90
## 3    2b       NA
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## 
## Parameters:
##         Estimate Std. Error t value Pr(>|t|)    
## tau     0.865508   0.243990   3.547 0.000393 ***
## alpha   0.877804   0.063130  13.905  < 2e-16 ***
## A     382.711524  46.064067   8.308  < 2e-16 ***
## k     181.086971  27.994519   6.469  1.1e-10 ***
## p      -0.021832   0.006205  -3.518 0.000439 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5791 on 4357 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 8.562e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   4357     1488.5                                
## 2   4356     1424.2  1 64.342   196.8 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2a 46690.90
## 2     4 46772.49
## 3     5 46581.74
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.00263    0.25225   3.975 7.16e-05 ***
## alpha   0.92051    0.06062  15.184  < 2e-16 ***
## a      12.37804    1.27551   9.704  < 2e-16 ***
## b     109.05738    6.07768  17.944  < 2e-16 ***
## c     113.94899    6.97245  16.343  < 2e-16 ***
## d       1.18472    0.05679  20.861  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5718 on 4356 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

plot residuals

## 
## ------
##  Shapiro-Wilk normality test
## 
## data:  stdres
## W = 0.96113, p-value < 2.2e-16
## 
## 
## ------
## 
##  Runs Test
## 
## data:  as.factor(run)
## Standard Normal = -31.622, p-value < 2.2e-16
## alternative hypothesis: two.sided

predict and plot

plotting 2

223 - Central Interior Broadleaf Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7425     1362.9                                
## 2   7424     1316.5  1 46.482  262.13 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 77578.06
## 2     2 77322.31
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.04382    0.14030    7.44 1.12e-13 ***
## alpha   0.79468    0.04635   17.15  < 2e-16 ***
## A     229.08813    9.44538   24.25  < 2e-16 ***
## k      88.54009    4.45866   19.86  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4211 on 7424 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 1.902e-06
##   (3 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq   Df  Sum Sq F value    Pr(>F)    
## 1   7424     1316.5                                   
## 2   7423     1312.4    1   4.078 23.0681 1.594e-06 ***
## 3   7423     1296.2    0   0.000                      
## 4   6380     1128.1 1043 168.074  0.9114    0.9733    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 77322.31
## 2    2a 77301.26
## 3    2b 77208.80
## 4    2c 66305.16
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.85759    0.12173   7.045 2.05e-12 ***
## alpha   0.71209    0.03381  21.062  < 2e-16 ***
## A     127.11413    4.61707  27.531  < 2e-16 ***
## k      37.02397    1.25609  29.475  < 2e-16 ***
## p       0.14094    0.01883   7.486 8.09e-14 ***
## s       2.17271    0.15279  14.220  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4205 on 6380 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 4.063e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7423     1331.0                                
## 2   7422     1284.9  1 46.068   266.1 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2c 66305.16
## 2     4 77405.77
## 3     5 77146.11
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.85759    0.12173   7.045 2.05e-12 ***
## alpha   0.71209    0.03381  21.062  < 2e-16 ***
## A     127.11413    4.61707  27.531  < 2e-16 ***
## k      37.02397    1.25609  29.475  < 2e-16 ***
## p       0.14094    0.01883   7.486 8.09e-14 ***
## s       2.17271    0.15279  14.220  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4205 on 6380 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 4.063e-06
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

231 - Southeastern Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   9262     2333.1                                
## 2   9261     2269.7  1 63.327  258.39 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 101389.5
## 2     2 101136.6
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.58614    0.16735   9.478   <2e-16 ***
## alpha   0.71207    0.04222  16.865   <2e-16 ***
## A     286.74176   11.66804  24.575   <2e-16 ***
## k      85.62977    3.37714  25.356   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4951 on 9261 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 1.425e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq   Df   Sum Sq  F value    Pr(>F)    
## 1   9261     2269.7                                     
## 2   9260     2263.6    1    6.179  25.2788 5.054e-07 ***
## 3   7771     2566.4 1489 -302.856  -0.6159         1    
## 4   7770     2496.6    1   69.841 217.3639 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model       AIC
## 1     2 101136.58
## 2    2a 101113.33
## 3    2b  85459.40
## 4    2c  85246.85
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau   1.506e+00  1.836e-01   8.201 2.76e-16 ***
## alpha 6.512e-01  2.258e-02  28.840  < 2e-16 ***
## A     1.450e+02  6.248e+00  23.205  < 2e-16 ***
## k     3.016e+01  1.222e+00  24.686  < 2e-16 ***
## p     1.076e-01  7.674e-03  14.024  < 2e-16 ***
## s     1.963e+00  9.744e-02  20.143  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5668 on 7770 degrees of freedom
## 
## Number of iterations to convergence: 12 
## Achieved convergence tolerance: 5.033e-06
##   (14 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   9260     2303.2                                
## 2   9259     2236.8  1 66.418  274.94 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model       AIC
## 1    2c  85246.85
## 2     4 101274.10
## 3     5 101004.99
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau   1.506e+00  1.836e-01   8.201 2.76e-16 ***
## alpha 6.512e-01  2.258e-02  28.840  < 2e-16 ***
## A     1.450e+02  6.248e+00  23.205  < 2e-16 ***
## k     3.016e+01  1.222e+00  24.686  < 2e-16 ***
## p     1.076e-01  7.674e-03  14.024  < 2e-16 ***
## s     1.963e+00  9.744e-02  20.143  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5668 on 7770 degrees of freedom
## 
## Number of iterations to convergence: 12 
## Achieved convergence tolerance: 5.033e-06
##   (14 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

232 - Outer Coastal Plain Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   9495     3706.9                                
## 2   9494     3585.7  1 121.15  320.78 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 106844.7
## 2     2 106531.1
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.12108    0.18182   6.166 7.29e-10 ***
## alpha   0.74475    0.03909  19.052  < 2e-16 ***
## A     303.18915   15.39059  19.700  < 2e-16 ***
## k      91.43033    4.53472  20.162  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6146 on 9494 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 6.341e-06
##   (8 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   9494     3585.7                                
## 2   9493     3580.8  1 4.8991  12.988 0.0003151 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model       AIC
## 1     2 106531.09
## 2    2a 106520.11
## 3    2b  88567.74
## 4    2c        NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.83554    0.18434   4.533 5.91e-06 ***
## alpha   0.58429    0.02782  21.001  < 2e-16 ***
## A     700.20538  245.84610   2.848  0.00441 ** 
## k     453.42710  284.16976   1.596  0.11061    
## s       0.77872    0.03365  23.145  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7037 on 7903 degrees of freedom
## 
## Number of iterations to convergence: 12 
## Achieved convergence tolerance: 2.783e-06
##   (32 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   9493     3686.3                                
## 2   9492     3559.9  1 126.38  336.99 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model       AIC
## 1    2b  88567.74
## 2     4 106795.73
## 3     5 106466.38
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.83554    0.18434   4.533 5.91e-06 ***
## alpha   0.58429    0.02782  21.001  < 2e-16 ***
## A     700.20538  245.84610   2.848  0.00441 ** 
## k     453.42710  284.16976   1.596  0.11061    
## s       0.77872    0.03365  23.145  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7037 on 7903 degrees of freedom
## 
## Number of iterations to convergence: 12 
## Achieved convergence tolerance: 2.783e-06
##   (32 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

234 - Lower Mississippi Riverine Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value   Pr(>F)    
## 1   1083     356.24                               
## 2   1082     339.64  1 16.601  52.886 6.77e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 12290.56
## 2     2 12240.74
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau      0.9025     0.6136   1.471   0.1417    
## alpha    0.8625     0.1098   7.854 9.68e-15 ***
## A     1162.2886   522.4216   2.225   0.0263 *  
## k      500.8135   237.4249   2.109   0.0351 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5603 on 1082 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 2.518e-06
##   (9 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in nls(get(paste("f_", Mod.Sel1, "a", sep = "")), data = G_234,  : 
##   step factor 0.000488281 reduced below 'minFactor' of 0.000976562
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq  Df Sum Sq F value    Pr(>F)    
## 1   1082     339.64                                 
## 2    823     192.48 259 147.16  2.4293 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model       AIC
## 1     2 12240.738
## 2    2a        NA
## 3    2b        NA
## 4    2c  8996.781
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.39536    0.41113   0.962    0.337    
## alpha   0.78925    0.07021  11.242  < 2e-16 ***
## A     395.47217  189.29921   2.089    0.037 *  
## k      99.96670   70.69953   1.414    0.158    
## p      -0.01220    0.01423  -0.857    0.392    
## s       1.13109    0.24438   4.628 4.28e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4836 on 823 degrees of freedom
## 
## Number of iterations to convergence: 40 
## Achieved convergence tolerance: 3.1e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   1081     355.31                                
## 2   1080     337.44  1 17.869  57.192 8.426e-14 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model       AIC
## 1    2c  8996.781
## 2     4 12291.729
## 3     5 12237.690
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.39536    0.41113   0.962    0.337    
## alpha   0.78925    0.07021  11.242  < 2e-16 ***
## A     395.47217  189.29921   2.089    0.037 *  
## k      99.96670   70.69953   1.414    0.158    
## p      -0.01220    0.01423  -0.857    0.392    
## s       1.13109    0.24438   4.628 4.28e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4836 on 823 degrees of freedom
## 
## Number of iterations to convergence: 40 
## Achieved convergence tolerance: 3.1e-06
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

plot residuals

## 
## ------
##  Shapiro-Wilk normality test
## 
## data:  stdres
## W = 0.95371, p-value = 1.756e-15
## 
## 
## ------
## 
##  Runs Test
## 
## data:  as.factor(run)
## Standard Normal = -10.597, p-value < 2.2e-16
## alternative hypothesis: two.sided

predict and plot

plotting 2

242 - Pacific Lowland Mixed Forest

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

plot residuals

## [1] "cannot plot residuals"

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

251 - Prairie Parkland (Temperate)

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   1884     530.49                                
## 2   1883     520.40  1 10.091  36.511 1.825e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 19835.95
## 2     2 19801.71
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.18712    0.25641   0.730    0.466    
## alpha   0.63418    0.09964   6.364 2.45e-10 ***
## A     298.73479   32.98617   9.056  < 2e-16 ***
## k     126.75170   16.66267   7.607 4.41e-14 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5257 on 1883 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 2.294e-06
##   (3 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df  Sum Sq F value Pr(>F)
## 1   1883     520.40                          
## 2   1882     520.18  1 0.21527  0.7788 0.3776
##   model      AIC
## 1     2 19801.71
## 2    2a 19802.93
## 3    2b 19791.60
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.2497     0.2619   0.954     0.34    
## alpha   0.6703     0.1001   6.699 2.76e-11 ***
## A     167.9953    18.6424   9.011  < 2e-16 ***
## k      47.7477     6.8077   7.014 3.22e-12 ***
## s       1.4586     0.1242  11.742  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5242 on 1882 degrees of freedom
## 
## Number of iterations to convergence: 9 
## Achieved convergence tolerance: 7.144e-06
##   (3 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   1882     514.77                                
## 2   1881     503.84  1 10.929    40.8 2.119e-10 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 19791.60
## 2     4 19783.19
## 3     5 19744.70
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.20201    0.25282   0.799    0.424    
## alpha   0.65588    0.09727   6.743 2.06e-11 ***
## a      21.61048    2.71796   7.951 3.16e-15 ***
## b     100.87424    7.44340  13.552  < 2e-16 ***
## c     107.78366   10.02239  10.754  < 2e-16 ***
## d       1.20475    0.09772  12.328  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5176 on 1881 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (3 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

255 - Prairie Parkland (Subtropical)

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    633     274.88                                
## 2    632     263.69  1 11.189  26.816 3.012e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 6532.531
## 2     2 6508.102
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.1250     0.5533   0.226 0.821390    
## alpha   0.7041     0.1238   5.687 1.98e-08 ***
## A     296.3420    70.9096   4.179 3.34e-05 ***
## k     128.7173    33.4109   3.853 0.000129 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6459 on 632 degrees of freedom
## 
## Number of iterations to convergence: 8 
## Achieved convergence tolerance: 3.502e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1    632     263.69                         
## 2    631     262.60  1 1.0917  2.6232 0.1058
##   model      AIC
## 1     2 6508.102
## 2    2a 6507.463
## 3    2b       NA
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     0.08586    0.54208   0.158  0.87420    
## alpha   0.73491    0.12324   5.963 4.11e-09 ***
## A     242.91208   59.23770   4.101 4.66e-05 ***
## k      88.31753   30.37767   2.907  0.00377 ** 
## p      -0.01872    0.01659  -1.128  0.25970    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6451 on 631 degrees of freedom
## 
## Number of iterations to convergence: 13 
## Achieved convergence tolerance: 6.185e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    631     262.85                                
## 2    630     247.58  1 15.265  38.844 8.404e-10 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2a 6507.463
## 2     4 6508.057
## 3     5 6472.004
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau   -0.08353    0.48167  -0.173    0.862    
## alpha  0.79872    0.11384   7.016 5.91e-12 ***
## a     11.01478    2.30217   4.785 2.14e-06 ***
## b     88.30672   10.55299   8.368 3.79e-16 ***
## c     58.54755    5.05845  11.574  < 2e-16 ***
## d      1.06975    0.09500  11.261  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6269 on 630 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

261 - California Coastal Chaparral Forest and Shrub

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit

  • add s model: does not fit

  • add s+p model: does not fit

  • unable to fit model (only 64 observations)

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

262 - California Dry Steppe

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit

  • add s model: does not fit

  • add s+p model: does not fit

  • unable to fit model (0 observations)

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

263 - California Coastal Steppe - Mixed Forest and Redwood Forest

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

313 - Colorado Plateau Semi-Desert

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

315 - Southwest Plateau and Plains Dry Steppe and Shrub

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

321 - Chihuahuan Semi-Desert

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

322 - American Semidesert and Desert

model selection 1

## Error in nls(f_1, data = G_322, start = c(tau = tau.start, A = A.start,  : 
##   step factor 0.000488281 reduced below 'minFactor' of 0.000976562
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
##   model AIC
## 1     1  NA
## 2     2  NA
## Warning in min(AIC1_322$AIC, na.rm = T): no non-missing arguments to min;
## returning Inf
## Error in h(simpleError(msg, call)) : 
##   error in evaluating the argument 'object' in selecting a method for function 'summary': object 'nls_322.' not found

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"
  • Cannot fit model
  • not enough data (only 3 observations)

331 - Great Plains/Palouse Dry Steppe

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"
  • Cannot fit model

332 - Great Plains Steppe

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value Pr(>F)
## 1    216     145.31                         
## 2    215     143.71  1 1.5947  2.3858 0.1239
##   model      AIC
## 1     1 2335.903
## 2     2 2335.487
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)  
## tau      0.5363     1.8691   0.287    0.774  
## alpha    0.5955     0.3508   1.698    0.091 .
## A     1637.2126  6225.7800   0.263    0.793  
## k     1408.2568  5450.4162   0.258    0.796  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8176 on 215 degrees of freedom
## 
## Number of iterations to convergence: 8 
## Achieved convergence tolerance: 2.297e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in nls(get(paste("f_", Mod.Sel1, "a", sep = "")), data = G_332,  : 
##   step factor 0.000488281 reduced below 'minFactor' of 0.000976562
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
##   Res.Df Res.Sum Sq Df Sum Sq F value  Pr(>F)  
## 1    215    143.711                            
## 2    147     83.108 68 60.603  1.5764 0.01166 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 2335.487
## 2    2a       NA
## 3    2b 1636.402
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)   
## tau     2.8328     3.9193   0.723  0.47097   
## alpha   0.9760     0.3039   3.211  0.00162 **
## A      85.6246    59.9298   1.429  0.15520   
## k      44.9763    23.3407   1.927  0.05592 . 
## s       1.7214     0.6592   2.611  0.00995 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7519 on 147 degrees of freedom
## 
## Number of iterations to convergence: 33 
## Achieved convergence tolerance: 8.579e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: does not fit
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value  Pr(>F)  
## 1    214     137.79                            
## 2    213     135.66  1 2.1317  3.3469 0.06873 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 1636.402
## 2     4 2328.276
## 3     5 2326.862
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)   
## tau     2.8328     3.9193   0.723  0.47097   
## alpha   0.9760     0.3039   3.211  0.00162 **
## A      85.6246    59.9298   1.429  0.15520   
## k      44.9763    23.3407   1.927  0.05592 . 
## s       1.7214     0.6592   2.611  0.00995 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7519 on 147 degrees of freedom
## 
## Number of iterations to convergence: 33 
## Achieved convergence tolerance: 8.579e-06
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

341 - Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit

  • log-normal phi model: does not fit

  • model not fitted because only 62 observations

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

342 - Intermountain Semi-Desert

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

411 - Everglades

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

M211 - Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   4834     775.91                                
## 2   4833     746.70  1  29.21  189.06 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 51220.60
## 2     2 51036.99
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.81073    0.28200   6.421 1.48e-10 ***
## alpha   0.87838    0.06026  14.577  < 2e-16 ***
## A     332.52443   23.77680  13.985  < 2e-16 ***
## k     173.44970   11.22937  15.446  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3931 on 4833 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 4.621e-06
##   (1 observation deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df Sum Sq F value   Pr(>F)    
## 1   4833     746.70                               
## 2   4832     738.81  1 7.8916  51.613 7.79e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 51036.99
## 2    2a 50987.60
## 3    2b 50965.66
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.52010    0.25515   5.958 2.74e-09 ***
## alpha   0.91818    0.06060  15.152  < 2e-16 ***
## A     181.37233   11.65322  15.564  < 2e-16 ***
## k      57.00052    4.15035  13.734  < 2e-16 ***
## s       1.46437    0.06059  24.167  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3901 on 4832 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 6.556e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   4832     765.27                                
## 2   4831     733.83  1 31.442  206.99 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 50965.66
## 2     4 51157.82
## 3     5 50956.90
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.56244    0.25921   6.028 1.79e-09 ***
## alpha   0.92196    0.06021  15.312  < 2e-16 ***
## a      12.85464    2.78039   4.623 3.87e-06 ***
## b     125.02105    8.62168  14.501  < 2e-16 ***
## c     166.06497   16.00981  10.373  < 2e-16 ***
## d       1.54602    0.09924  15.579  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3897 on 4831 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

plot residuals

## 
## ------
##  Shapiro-Wilk normality test
## 
## data:  stdres
## W = 0.98658, p-value < 2.2e-16
## 
## 
## ------
## 
##  Runs Test
## 
## data:  as.factor(run)
## Standard Normal = -33.203, p-value < 2.2e-16
## alternative hypothesis: two.sided

predict and plot

plotting 2

M221 - Eastern Broadleaf Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7196     1178.4                                
## 2   7195     1129.9  1 48.537  309.08 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 79716.90
## 2     2 79416.11
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.9250     0.1344   6.881 6.43e-12 ***
## alpha   0.9052     0.0490  18.474  < 2e-16 ***
## A     258.3782     8.9769  28.783  < 2e-16 ***
## k      63.0567     2.7132  23.241  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3963 on 7195 degrees of freedom
## 
## Number of iterations to convergence: 7 
## Achieved convergence tolerance: 4.514e-06
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7195     1129.9                                
## 2   7194     1124.6  1 5.3301  34.098 5.469e-09 ***
## 3   7194     1114.1  0 0.0000                      
## 4   7193     1104.5  1 9.6107  62.590 2.928e-15 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 79416.11
## 2    2a 79384.07
## 3    2b 79316.88
## 4    2c 79256.50
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.08137    0.14102   7.668 1.98e-14 ***
## alpha   0.91493    0.04846  18.879  < 2e-16 ***
## A     159.40185    5.13905  31.018  < 2e-16 ***
## k      35.98938    0.89911  40.028  < 2e-16 ***
## p       0.16407    0.01612  10.178  < 2e-16 ***
## s       2.45792    0.14847  16.555  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3919 on 7193 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 9.267e-06
##   (2 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   7194     1150.4                                
## 2   7193     1101.1  1 49.273  321.88 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2c 79256.50
## 2     4 79547.40
## 3     5 79234.25
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau     1.09590    0.14141    7.75 1.05e-14 ***
## alpha   0.91379    0.04844   18.86  < 2e-16 ***
## a      25.29252    2.30219   10.99  < 2e-16 ***
## b     123.73765    4.27083   28.97  < 2e-16 ***
## c     104.02553    3.50472   29.68  < 2e-16 ***
## d       1.33508    0.05096   26.20  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3913 on 7193 degrees of freedom
## 
## Algorithm "port", convergence message: relative convergence (4)
##   (2 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M223 - Ozark Broadleaf Forest Meadow

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    750     92.017                                
## 2    749     84.644  1 7.3729  65.242 2.646e-15 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 7647.094
## 2     2 7586.205
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.1367     0.2902   0.471    0.638    
## alpha   0.9558     0.1105   8.651  < 2e-16 ***
## A     284.8395    36.3007   7.847 1.48e-14 ***
## k      93.3971    16.0340   5.825 8.48e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3362 on 749 degrees of freedom
## 
## Number of iterations to convergence: 5 
## Achieved convergence tolerance: 3.457e-07
##   (2 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in nls(get(paste("f_", Mod.Sel1, "c", sep = "")), data = G_M223,  : 
##   number of iterations exceeded maximum of 50
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
##   Res.Df Res.Sum Sq Df  Sum Sq F value  Pr(>F)  
## 1    749     84.644                             
## 2    748     84.115  1 0.52843  4.6991 0.03049 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 7586.205
## 2    2a 7583.489
## 3    2b 6048.497
## 4    2c       NA
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau    -0.02502    0.27448  -0.091  0.92739    
## alpha   0.90237    0.09777   9.230  < 2e-16 ***
## A     214.41085   57.48780   3.730  0.00021 ***
## k      53.39305   23.09154   2.312  0.02111 *  
## s       1.22219    0.24916   4.905 1.21e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3477 on 594 degrees of freedom
## 
## Number of iterations to convergence: 9 
## Achieved convergence tolerance: 6.551e-06
##   (3 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: does not fit

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    748     91.700                                
## 2    747     84.333  1 7.3666  65.251 2.643e-15 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 6048.497
## 2     4 7648.498
## 3     5 7587.438
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##        Estimate Std. Error t value Pr(>|t|)    
## tau    -0.02502    0.27448  -0.091  0.92739    
## alpha   0.90237    0.09777   9.230  < 2e-16 ***
## A     214.41085   57.48780   3.730  0.00021 ***
## k      53.39305   23.09154   2.312  0.02111 *  
## s       1.22219    0.24916   4.905 1.21e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.3477 on 594 degrees of freedom
## 
## Number of iterations to convergence: 9 
## Achieved convergence tolerance: 6.551e-06
##   (3 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M231 - Ouachita Mixed Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    836     173.53                                
## 2    835     170.37  1 3.1556  15.466 9.099e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 8690.011
## 2     2 8676.613
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     1.0845     0.6144   1.765   0.0779 .  
## alpha   0.6828     0.1673   4.082 4.89e-05 ***
## A     336.2851    72.6004   4.632 4.20e-06 ***
## k     183.2880    41.2191   4.447 9.90e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4517 on 835 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 8.727e-06
##   (1 observation deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq  Df  Sum Sq F value Pr(>F)
## 1    835     170.37                           
## 2    834     170.22   1  0.1486  0.7280 0.3938
## 3    674     142.26 160 27.9596  0.8279 0.9282
## 4    673     142.26   1  0.0069  0.0327 0.8564
##   model      AIC
## 1     2 8676.613
## 2    2a 8677.881
## 3    2b 6981.181
## 4    2c 6983.148
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.5623     0.5070   1.109    0.268    
## alpha   0.8573     0.1138   7.531 1.63e-13 ***
## A     347.4775   233.5990   1.487    0.137    
## k     173.4408   187.9875   0.923    0.357    
## s       0.9738     0.1591   6.122 1.57e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4594 on 674 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 1.693e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    834     173.48                                
## 2    833     170.22  1 3.2592  15.949 7.081e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1    2b 6981.181
## 2     4 8693.783
## 3     5 8679.871
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     0.5623     0.5070   1.109    0.268    
## alpha   0.8573     0.1138   7.531 1.63e-13 ***
## A     347.4775   233.5990   1.487    0.137    
## k     173.4408   187.9875   0.923    0.357    
## s       0.9738     0.1591   6.122 1.57e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.4594 on 674 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 1.693e-06
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M242 - Cascade Mixed Forest

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

## [1] "cannot plot observed vs. predicted"

M261 - Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   1870     1216.8                                
## 2   1869     1208.6  1 8.2003  12.681 0.0003786 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 23583.09
## 2     2 23572.43
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     5.6743     3.1629   1.794 0.072978 .  
## alpha   0.4630     0.1238   3.741 0.000189 ***
## A     242.4795    82.5212   2.938 0.003340 ** 
## k     231.8338    32.1288   7.216 7.76e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.8041 on 1869 degrees of freedom
## 
## Number of iterations to convergence: 6 
## Achieved convergence tolerance: 5.895e-06
##   (6 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Error in numericDeriv(form[[3L]], names(ind), env, central = nDcentral) : 
##   Missing value or an infinity produced when evaluating the model
## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq   Df  Sum Sq F value    Pr(>F)    
## 1   1869    1208.57                                   
## 2   1868    1192.61    1   15.96 25.0048 6.254e-07 ***
## 3    324      93.71 1544 1098.90  2.4607 < 2.2e-16 ***
## 4    323      93.31    1    0.40  1.3989    0.2378    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model       AIC
## 1     2 23572.426
## 2    2a 23549.521
## 3    2b  3485.785
## 4    2c  3486.363
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     1.8945     1.2923   1.466   0.1436    
## alpha   0.6580     0.2763   2.381   0.0178 *  
## A      65.9248    13.3842   4.926 1.34e-06 ***
## k      30.8188     3.7464   8.226 4.75e-15 ***
## s      42.5340   187.7406   0.227   0.8209    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5378 on 324 degrees of freedom
## 
## Number of iterations to convergence: 21 
## Achieved convergence tolerance: 8.771e-06
##   (1 observation deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1   1868     1197.2                                
## 2   1867     1188.6  1 8.6142  13.531 0.0002413 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model       AIC
## 1    2b  3485.785
## 2     4 23556.736
## 3     5 23545.211
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     1.8945     1.2923   1.466   0.1436    
## alpha   0.6580     0.2763   2.381   0.0178 *  
## A      65.9248    13.3842   4.926 1.34e-06 ***
## k      30.8188     3.7464   8.226 4.75e-15 ***
## s      42.5340   187.7406   0.227   0.8209    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5378 on 324 degrees of freedom
## 
## Number of iterations to convergence: 21 
## Achieved convergence tolerance: 8.771e-06
##   (1 observation deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M262 - California coastal range - coniferous forest - open woodland - shrub meadow

Model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit

  • log-normal phi model: does not fit

  • model can fit - but K is negative (only 19 observations) - model excluded

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M313 - Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M331 - Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M332 - Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M333 - Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2

M334 - Black Hills Coniferous Forest

model selection 1

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    308     124.65                                
## 2    307     115.58  1 9.0627  24.071 1.509e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     1 3119.284
## 2     2 3097.808
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     1.6153     1.3279   1.216  0.22476    
## alpha   0.8352     0.1541   5.422 1.20e-07 ***
## A      69.6423    17.1618   4.058 6.28e-05 ***
## k      44.9466    15.3373   2.931  0.00364 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6136 on 307 degrees of freedom
## 
## Number of iterations to convergence: 4 
## Achieved convergence tolerance: 9.238e-06
##   (5 observations deleted due to missingness)

summary

  • simple model: fits
  • alpha model: fits

model selection 2

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2/(k + STDAGE_t2)))
## Model 3: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (A * STDAGE_t2^s/(k^s + STDAGE_t2^s))
## Model 4: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (p * A + ((1 - p) * A * STDAGE_t2^s/(k^s + STDAGE_t2^s)))
##   Res.Df Res.Sum Sq Df  Sum Sq F value Pr(>F)
## 1    307     115.58                          
## 2    306     115.48  1 0.10112  0.2679 0.6051
## 3    306     114.90  0 0.00000               
## 4    305     114.63  1 0.26693  0.7102 0.4000
##   model      AIC
## 1     2 3097.808
## 2    2a 3099.536
## 3    2b 3097.962
## 4    2c 3099.239
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     1.6153     1.3279   1.216  0.22476    
## alpha   0.8352     0.1541   5.422 1.20e-07 ***
## A      69.6423    17.1618   4.058 6.28e-05 ***
## k      44.9466    15.3373   2.931  0.00364 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6136 on 307 degrees of freedom
## 
## Number of iterations to convergence: 4 
## Achieved convergence tolerance: 9.238e-06
##   (5 observations deleted due to missingness)

summary

  • add p model: fits
  • add s model: fits
  • add s+p model: fits

model selection 3

## Analysis of Variance Table
## 
## Model 1: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
## Model 2: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * B_L_prop) * (a + b * exp(-((log(STDAGE_t2/c))/d)^2))
##   Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
## 1    306     123.46                                
## 2    305     115.54  1 7.9259  20.923 6.962e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   model      AIC
## 1     2 3097.808
## 2     4 3120.315
## 3     5 3101.680
## 
## Formula: B_plt_t2_MgHa ~ (1 + (MEASTIME_t2 - 1990) * tau/100) * (1 - alpha * 
##     B_L_prop) * A * STDAGE_t2/(k + STDAGE_t2)
## 
## Parameters:
##       Estimate Std. Error t value Pr(>|t|)    
## tau     1.6153     1.3279   1.216  0.22476    
## alpha   0.8352     0.1541   5.422 1.20e-07 ***
## A      69.6423    17.1618   4.058 6.28e-05 ***
## k      44.9466    15.3373   2.931  0.00364 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6136 on 307 degrees of freedom
## 
## Number of iterations to convergence: 4 
## Achieved convergence tolerance: 9.238e-06
##   (5 observations deleted due to missingness)

summary

  • simple log-normal model: fits
  • log-normal phi model: fits

predict and plot

plotting 2

M341 - Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow

model selection 1

summary

  • simple model: does not fit
  • alpha model: does not fit

model selection 2

summary

  • add p model: does not fit
  • add s model: does not fit
  • add s+p model: does not fit

model selection 3

summary

  • simple log-normal model: does not fit
  • log-normal phi model: does not fit

predict and plot

## [1] "cannot plot data with prediction"

plotting 2


Fitted parameters

Best / selected models by ecoprovince

Code Ecoregion Sel.Mod.2 Sel.Mod.3 Best.Mod
211 Northeastern Mixed Forest 2b 5 5
212 Laurentian Mixed Forest 2b 5 5
221 Eastern Broadleaf Forest 2b 2b 2b
222 Midwest Broadleaf Forest 2a 5 5
223 Central Interior Broadleaf Forest 2c 2c 2c
231 Southeastern Mixed Forest 2c 2c 2c
232 Outer Coastal Plain Mixed Forest 2b 2b 2b
234 Lower Mississippi Riverine Forest 2c 2c 2c
242 Pacific Lowland Mixed Forest NA NA NA
251 Prairie Parkland (Temperate) 2b 5 5
255 Prairie Parkland (Subtropical) 2a 5 5
261 California Coastal Chaparral Forest and Shrub NA NA NA
262 California Dry Steppe NA NA NA
263 California Coastal Steppe - Mixed Forest and Redwood Forest NA NA NA
313 Colorado Plateau Semi-Desert NA NA NA
315 Southwest Plateau and Plains Dry Steppe and Shrub NA NA NA
321 Chihuahuan Semi-Desert NA NA NA
322 American Semidesert and Desert NA NA NA
331 Great Plains/Palouse Dry Steppe NA NA NA
332 Great Plains Steppe 2b 2b 2b
341 Intermountain Semi-Desert and Desert NA NA NA
342 Intermountain Semi-Desert NA NA NA
411 Everglades NA NA NA
M211 Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow 2b 5 5
M221 Central Appalachian Broadleaf Forest - Coniferous Forest - Meadow 2c 5 5
M223 Ozark Broadleaf Forest Meadow 2b 2b 2b
M231 Ouachita Mixed Forest 2b 2b 2b
M242 Cascade Mixed Forest NA NA NA
M261 Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow 2b 2b 2b
M262 California Coastal Range Coniferous Forest - Open Woodland - Shrub - Meadow NA NA NA
M313 Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow NA NA NA
M331 Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow NA NA NA
M332 Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow NA NA NA
M333 Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow NA NA NA
M334 Black Hills Coniferous Forest 2 2 2
M341 Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow NA NA NA

table by ecoprovince

Code Ecoregion region n.obs n.plots tau tau.variance tau.2.5 tau.97.5 alpha alpha.variance alpha.2.5 alpha.97.5 A A.2.5 A.97.5 k k.2.5 k.97.5 a a.2.5 a.97.5 b b.se b.2.5 b.97.5 c c.2.5 c.97.5 d d.2.5 d.97.5
211 Northeastern Mixed Forest east 4838 2419 1.5826572 0.0859229 1.0079948 2.1573197 1.1088849 0.0041371 0.9827878 1.2349819 214.23259 176.03157 252.43361 72.38593 55.136299 89.63556 23.88706 19.889474 27.88464 108.49559 NA 94.71761 122.27358 132.60613 115.29891 149.91335 1.194045 1.060323 1.327767
212 Laurentian Mixed Forest east 12962 6481 1.2146588 0.0184839 0.9481722 1.4811455 0.8239244 0.0013527 0.7518336 0.8960153 112.92651 104.53007 121.32295 41.69107 37.451306 45.93084 13.98018 12.400597 15.55976 71.24841 NA 67.14756 75.34926 116.92499 108.58279 125.26720 1.411748 1.327560 1.495935
221 Eastern Broadleaf Forest east 5446 2723 0.2338300 0.0169227 -0.0211930 0.4888529 0.8162353 0.0009875 0.7546312 0.8778394 293.53650 251.49470 335.57830 64.20469 51.624025 76.78536 NA NA NA NA NA NA NA NA NA NA NA NA NA
222 Midwest Broadleaf Forest east 3552 1776 1.0026303 0.0636318 0.5080853 1.4971753 0.9205052 0.0036751 0.8016535 1.0393569 382.71152 292.40252 473.02052 181.08697 126.203476 235.97047 12.37804 9.877379 14.87869 109.05738 NA 97.14204 120.97273 113.94899 100.27945 127.61854 1.184719 1.073380 1.296058
223 Central Interior Broadleaf Forest east 6388 3194 0.8575908 0.0148186 0.6189556 1.0962260 0.7120924 0.0011431 0.6458139 0.7783709 127.11413 118.06312 136.16514 37.02397 34.561601 39.48633 NA NA NA NA NA NA NA NA NA NA NA NA NA
231 Southeastern Mixed Forest east 7790 3895 1.5059413 0.0337173 1.1459915 1.8658910 0.6511538 0.0005098 0.6068943 0.6954133 144.97930 132.73182 157.22678 30.16077 27.765806 32.55574 NA NA NA NA NA NA NA NA NA NA NA NA NA
232 Outer Coastal Plain Mixed Forest east 7940 3970 0.8355407 0.0339820 0.4741816 1.1968998 0.5842890 0.0007741 0.5297507 0.6388272 700.20538 218.28207 1182.12870 453.42710 -103.620712 1010.47492 NA NA NA NA NA NA NA NA NA NA NA NA NA
234 Lower Mississippi Riverine Forest east 830 415 0.3953616 0.1690247 -0.4116176 1.2023407 0.7892546 0.0049287 0.6514526 0.9270566 395.47217 23.90610 767.03825 99.96670 -38.805927 238.73932 NA NA NA NA NA NA NA NA NA NA NA NA NA
242 Pacific Lowland Mixed Forest west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
251 Prairie Parkland (Temperate) east 1392 696 0.2020061 0.0639157 -0.2938224 0.6978346 0.6558780 0.0094616 0.4651081 0.8466478 167.99527 131.43342 204.55712 47.74774 34.396282 61.09920 21.61048 16.279939 26.94101 100.87424 NA 86.27605 115.47242 107.78366 88.12748 127.43984 1.204748 1.013088 1.396409
255 Prairie Parkland (Subtropical) east 444 222 -0.0835282 0.2320018 -1.0293926 0.8623362 0.7987237 0.0129602 0.5751659 1.0222815 242.91208 126.58518 359.23897 88.31753 28.663971 147.97109 11.01478 6.493934 15.53563 88.30672 NA 67.58342 109.03002 58.54755 48.61408 68.48103 1.069750 0.883204 1.256296
261 California Coastal Chaparral Forest and Shrub west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
262 California Dry Steppe west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
263 California Coastal Steppe - Mixed Forest and Redwood Forest west 4 2 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
313 Colorado Plateau Semi-Desert west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
315 Southwest Plateau and Plains Dry Steppe and Shrub west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
321 Chihuahuan Semi-Desert west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
322 American Semidesert and Desert west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
331 Great Plains/Palouse Dry Steppe west 118 59 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
332 Great Plains Steppe west 154 77 2.8327582 15.3607851 -4.9126609 10.5781773 0.9759876 0.0923741 0.3753487 1.5766264 85.62460 -32.81058 204.05978 44.97633 -1.150377 91.10303 NA NA NA NA NA NA NA NA NA NA NA NA NA
341 Intermountain Semi-Desert and Desert west 4 2 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
342 Intermountain Semi-Desert west 2 1 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
411 Everglades east 66 33 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M211 Adirondack-New England Mixed forest - Coniferous Forest - Alpine Meadow east 5108 2554 1.5624417 0.0671905 1.0542697 2.0706137 0.9219623 0.0036253 0.8039230 1.0400017 181.37233 158.52671 204.21795 57.00052 48.863953 65.13709 12.85464 7.403809 18.30547 125.02105 NA 108.11864 141.92346 166.06497 134.67846 197.45149 1.546021 1.351465 1.740576
M221 Central Appalachian Broadleaf Forest - Coniferous Forest - Meadow east 5186 2593 1.0959008 0.0199964 0.8186986 1.3731030 0.9137907 0.0023469 0.8188253 1.0087562 159.40185 149.32781 169.47590 35.98938 34.226867 37.75190 25.29252 20.779542 29.80549 123.73765 NA 115.36557 132.10973 104.02553 97.15524 110.89581 1.335080 1.235176 1.434983
M223 Ozark Broadleaf Forest Meadow east 602 301 -0.0250231 0.0753381 -0.5640881 0.5140420 0.9023719 0.0095586 0.7103585 1.0943853 214.41085 101.50679 327.31492 53.39305 8.042058 98.74405 NA NA NA NA NA NA NA NA NA NA NA NA NA
M231 Ouachita Mixed Forest east 680 340 0.5622631 0.2570724 -0.4332717 1.5577979 0.8572806 0.0129590 0.6337615 1.0807997 347.47751 -111.19185 806.14687 173.44078 -195.670788 542.55234 NA NA NA NA NA NA NA NA NA NA NA NA NA
M242 Cascade Mixed Forest west 34 17 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M261 Sierran Steppe - Mixed Forest - Coniferous Forest - Alpine Meadow west 330 165 1.8945480 1.6701111 -0.6478649 4.4369609 0.6580466 0.0763610 0.1144095 1.2016837 65.92482 39.59387 92.25577 30.81882 23.448443 38.18920 NA NA NA NA NA NA NA NA NA NA NA NA NA
M262 California Coastal Range Coniferous Forest - Open Woodland - Shrub - Meadow west 8 4 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M313 Arizona-New Mexico Mountains Semi-Desert - Open Woodland - Coniferous Forest - Alpine Meadow west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M331 Southern Rocky Mountain Steppe - Open Woodland - Coniferous Forest - Alpine Meadow west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M332 Middle Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow west 20 10 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M333 Northern Rocky Mountain Steppe - Coniferous Forest - Alpine Meadow west 22 11 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
M334 Black Hills Coniferous Forest west 306 153 1.6153148 1.7633679 -0.9976589 4.2282885 0.8351977 0.0237321 0.5320656 1.1383298 69.64228 35.87271 103.41185 44.94656 14.767046 75.12608 NA NA NA NA NA NA NA NA NA NA NA NA NA
M341 Nevada-Utah Mountains Semi-Desert - Coniferous Forest - Alpine Meadow west 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA

plot tau

map

## OGR data source with driver: ESRI Shapefile 
## Source: "C:\Users\hogan.jaaron\Dropbox\FIA_R\Mapping\S_USA.EcoMapProvinces\S_USA.EcoMapProvinces.shp", layer: "S_USA.EcoMapProvinces"
## with 37 features
## It has 17 fields
## Integer64 fields read as strings:  PROVINCE_ PROVINCE_I

plot alpha (biomass compensation effect)

plot A (asymptote of B)

## Warning: Removed 19 rows containing missing values (`geom_point()`).

plot k (stand age at half biomass asymptote)

## Warning: Removed 19 rows containing missing values (`geom_point()`).